Superfluidity of indirect excitons in a quantum dot. 
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The superfluidity and Bose-Einstein condensation of indirect excitons in 
two-dimensional quantum dot are studied by path-integral Monte Carlo simu- 
lations. The temperature dependence of superfluid and Bose-condensed frac- 
tion are calculated at different strengths of interaction. Using the Kosterlitz- 
Thouless recursion relations, we also predict behavior of superfluid fraction 
in macroscopic large systems. 
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I. INTRODUCTION 

Low- dimensional electron-hole systems in semiconductor heterostructures, quantum dots 
and quantum wires are of increasing experimental and theoretical attention now. Concep- 
tually new physical effects can take place in these systems.^'! One of the most interesting 
examples of systems considered is the structure that consists of two vertically connected 
two-dimensional quantum wells. Under the action of laser radiation, each well is occupied 
by carriers of different types, which form a two-layered electron-hole system. As has been 
shown in Ref. |3], superfluidity of indirect excitons, which are bounded states of pairs of 
spatially-separated electrons and holes, can appear in the system and manifests itself as 
persistent electric currents in each two-dimensional layer. The behavior of phase diagram in 
strong; magnetic field! and Josephson effects! have been also studied. A number of interest- 
ing experiments for magnetoexcitons with spatially-separated carriers in coupled quantum 
wells have been recently carried outJUl 

The surface roughness can lead to the excitons localization^ and thus it behaves like a 
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natural quantum dot.i Excitons can also be confined in vertically coupled artificial quanti 
dots.0 Hence, the investigation of "lakes" of excitons is essential. 

The goal of this article is to study the onset of the superfluidity of indirect excitons in 
two-dimensional (2D) quantum dot. In a 2D system at low densities n <C a~ 2 (a is exciton 
radius), potential barrier that is generated by dipole-dipole repulsion of indirect excitons 
suppress exchanges of constituent particles, electrons and holes.i Therefore the system of 
2D indirect excitons at low density can be treated as Bose system. (Other situations in which 
excitons can be assumed to be bosons were discussed in Refs. p|-|l3|.) Moreover, under the 
same conditions indirect excitons can be thought of as point-like particles with unidirected 
dipole momenta. Another system which can be mapped onto this model is that of trapped 
Bose atoms with dipole momenta induced by d.c. electric field or standing electromagnetic 
wave in resonance with atomic transition. For the sake of definiteness we present explanation 
below in terms of excitonic cluster. In the absence of interaction between particles, we use 
a semi-analytical method based on a recursion relation for canonical partition functions. 
Clusters of interacting excitons are explored by the path-integral Monte Carlo technique. 

The paper is organized as follows. In the next section the model and controlling dimen- 
sionless parameters are introduced. In Sec. |TTTj we outline the Monte Carlo technique and a 



semi-analytical approach. Sec. IV is devoted to the discussion of the results. 



II. THE MODEL HAMILTONIAN 



As discussed above, we treat indirect exciton as a structureless boson with effective mass 
m and effective dipole moment d perpendicular to the plane of 2D quantum dot. The 
confining potential of the quantum dot we take in a harmonic form. Hence, the model 
Hamiltonian is given by 
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Choosing as length unit tq = yfi/muj, we introduce dimensionless controlling parameters for 
the system at the finite temperature f3~ l : dimensionless temperature T = / {h 2 / (2mrQ)) 
and interaction constant k = (|d| 2 / Vq) / (h 2 / '{2mr'fy) . 

III. SIMULATION METHODS AND CALCULATED PROPERTIES 

The canonical ensemble of noninteracting Bose system (k = 0) can be studied with the 
help of simple semi-analytical procedure analogous to that has been used for particles in a 
harmonic trap interacted by a harmonic law.0 The method that has been used is based on 
universal connections between partition functions of canonical and grand canonical ensem- 
bles 
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Here \x is the chemical potential of the grand canonical ensemble and N is the number of 
particles in the canonical ensembles. Starting from Eq. (|2|), one can derive a recursion 
relation for the canonical partition functions of a noninteracting Bose gas with an arbitrary 
one-particle spectrum 
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Z n(P) = ^ E Zn(P)Zi((N - n)(3), (3) 
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Hence, the numerical procedure consists in calculations of the canonical one-particle partition 
functions at temperatures T/N, 2T/N, ... T and iterations of Eq. (|3|). Generalization 
of the Eq. (|3|) to the case of partition function derivatives, in form of which interesting 
thermodynamic averages (see this section below) can be presented, is obvious. 

Cluster of interacting excitons has been studied by the path-integral Monte Carlo cal- 
culations (for comprehensive review of application the path-integral method to the Bose 



systems see Ref. [L6D. This technique is based on the Trotter procedure as a result of which 
thermodynamic averages on D dimensional quantum system are expressed via averages on 
D + 1 dimensional classical system, which is created by multiplication of the configuration 
space of the initial system along imaginary time direction. Calculated in so called primitive 
approximation, the potential energy of the effective classical system that corresponds to the 
quantum model with Hamiltonian ([[]) is 
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where P is the number of effective system layers, N is the number of particles in one layer, 
which is equal to the number of particles in the original quantum system and r/^ stands 
for the 2D radius vector of i-th particle in the layer p. Bose statistics is taken into account 
by summation over particles permutations V during Monte Carlo simulation. @ The energy 
of the effective system depends on V by assumption r/ P ^ = r^s. 

Our attention has been focused on superfluid properties of excitons. One of the most 
dramatic manifestations of the superfluidity is the deviation of the effective moment of inertia 
from its rigid or "classical" value l c \ = (j^iLi mrf^ . By definition, the effective moment of 
inertia is the linear response function to a uniform rotation field with angular frequency Q 
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Only the normal fluid gives rise to the / e //-^ Thus, the superfluid fraction of the cluster is 
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In order to calculate v s in case of ideal gas in a trap (i.e. at k = 0, see definition of k 
near Eq. (0)) by using above mentioned recurrent procedure, firstly, one should express J e // 
and I c i in form of partition function derivatives 
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Here Z Q = tr |exp(— (3(H — f2L z )) j is the partition function of the system rotating with the 
frequency Q, Z = Z Q=0 is the partition function of the system at rest. Secondly, explicit 
expression for one-particle partition function Z n is needed. The spectrum of the Hamil- 
tonian H n = H — QL Z at k = is e n ,m = hu(2n + \m\ + 1) — hQn, where n — 0, 1, . . . and 
m = 0, ±1, ±2, . . . are radial and angular quantum numbers (hm is the eigenvalue of the 
operator L z ). After some algebra, one can write 
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Calculation of superfluid fraction in case of k ^ have been performed by path-integral 
Monte Carlo simulations. Applying the Trotter procedure to the quantum-statistics averages 
by which I d and I e ff are defined, one takesEH 
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Here averagings are performed on the effective classical system, the moments l c \ and J e // 
are measured in units of mrQ — h/u. 

The fraction of Bose-condensed particles we calculate by the method introduced in Ref. 



18: the maximum length of permutation cycle in path-integral Monte Carlo simulation which 



has a nonzero probability is the number of particles in the condensate. 



IV. DISCUSSION OF THE SIMULATION RESULTS 

Calculated temperature dependencies of superfluid fraction are shown in Fig. 1. 
One can see superfluidity suppression induced by interaction at a fixed dot angular 
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frequency. Condensate depletion is obvious too (see Fig. 2). However after defin- 
ing "average density" n = N/{2n < r 2 >) and plotting v s vs. dimensionless temperature 
a = m/(ph 2 n) = 2nI c i/(ph 2 N 2 ), i.e. at a fixed density, one can see (Fig. 3) some rising of 
the superffuid fraction with interaction. But interaction, of course, depletes condensate even 
at fixed density (Fig. 4). 

The investigated system is fully analogous to the atoms in a plane harmonic trap. BEC 
occurs for noninteracting system in two-dimensional harmonic trap in thermodynamic limit 
(contrary to the situation in homogeneous systems). For systems with repulsive interaction 
there is only some critical temperatures of transition to the new state without condensate 
(see Ref. [19] and references therein). 

Assuming Kosterlitz-Thouless (KT) scenario, we have estimated the critical temperature 
of the phase transition. Nonanalytic behavior of the free energy of the system and conse- 
quently a phase transition can be observed only in macroscopic systems. But in a computer 
simulation one deals with finite systems. Hence, the extrapolating procedure is necessary. 
We define thermodynamic limit as u — > 0, N — > oo at n = N/(2tc < r 2 >) = constant. In 
case of weak interaction this definition coincides with requirement Nuo 2 = constant. 

The extrapolation procedure (see Ref. 120) is based on mapping of a 2D system of vortices 



to the 2D logarithmic gas, which has KT transition.0H 

The temperature of this gas, which is expressed in units of square of logarithmic 
gas charge, is T cg = mj (2nh 2 n(3) = a/(2ir), the chemical potential (in the same units) is 
\i = —T cg j3E c , where E c is the vortex core energy in superfluid, and dielectric constant is 
e = \jv s . Kosterlitz-Thouless recursion relations lead to the universal jump of 1/ {T c ge) from 
4 to at the transition temperature. But in the finite systems superfluid vanishes smoothly 
enough (Figs. 1, 3). Such type of behavior can be analytically accounted by integrating 
recursion relations up to size of system (not to infinity) .0 

The parameters of curve u s (T cg ) are vortex core energy and diameter, which are defined by 
fitting analytic approximation to the Monte Carlo results. Then full integration of Kosterlitz- 
Thouless relations at fixed vortex core energy and core diameter gives critical value of 
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temperature. 

Scaling results for different values of interaction constant k gives, e.g., a^r(O.l) 
a KT {l) = 1.19. 



= 0.92, 
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FIG. I. Superfluid fraction v s vs. dimensionless temperature T = 2/(/3huj) for excitonic 
cluster with number of particles N = 37: solid line - k — (noninteracting case), solid 
triangles - k = 0.01, open circles - k = 0.1, solid circles - k — 1, open squares k = 2.38, 
solid squares k = 23.5, open triangle k = 48.25. 

FIG. 2. Condensed fraction u vs. T = 2/(/3hu) for excitonic cluster with number of 
particles iV = 37: solid line - k — (noninteracting case), open circles - k = 0.1, solid 
circles - k = 1. 

FIG. 3. Superfluid fraction v s vs. a = 2i[I c i/(f31i 2 N 2 : solid line - k — (noninteracting 
case), solid triangles - k = 0.01, open circles - k = 0.1, solid circles - k — 1, open squares - 
k = 2.38, solid squares - k — 23.5. 

FIG. 4. Condensed fraction u vs. a = 27rI d /(/3h 2 N 2 : solid line - k — (noninteracting 
case), open circles — A: = 0.1, solid circles - k — 1. 
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